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ABSTRACT 

We present the first results of a serendipitous search for clusters of galaxies in deep 
ROSAT-PSPC pointed observations at high galactic latitude. The survey is being carried 
out using a Wavelet based Detection Algorithm which is not biased against extended, low 
surface brightness sources. A new fiux-diameter limited sample of 10 cluster candidates 
has been created from ~ 3 deg^ surveyed area. Preliminary CCD observations have 
revealed that a large fraction of these candidates correspond to a visible enhancement in 
the galaxy surface density, and several others have been identified from other surveys. 
We believe these sources to be either low-moderate redshift groups or intermediate to 
high redshift clusters. We show X-ray and optical images of some of the clusters identified 
to date. We present, for the first time, the derived number density of the galaxy clusters 
to a fiux limit of 1 • 10"^'* erg cm~^ s~^ (0.5-2.0 keV). This extends the log A^-log of 
previous cluster surveys by more than one decade in fiux. Results are compared to 
theoretical predictions for cluster number counts. 

Subject headings: galaxies: clusters: general — X-rays: galaxies, general — techniques: 
image processing 
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1. Introduction 

Clusters of galaxies are powerful tools for study- 
ing the large scale structure and evolution of the uni- 
verse. Current data provide conflicting constraints 
on cluster evolution as well as on possible cosmologi- 
cal models depending on whether cluster samples are 
optically or X-ray selected (see e.g. the discussion 
in Briel & Henry, 1993). These discrepancies could 
reflect the well known problems which characterize 
optical studies, especially at high redshifts, such as 
chance coincidences of unvirialized systems, biases to- 
wards increasingly rich systems, and biases towards 
systems with anomalously bright galaxies. Optically 
selected X-ray observed cluster samples (Castander 
et al. 1994, Bower et al. 1994) still suffer from the 
biases of the original selection process. By using X- 
ray surveys (Gioia et al. 1990, Piccinotti et al. 1982) 
which have the advantage of selecting physical ob- 
jects, deep virialized potential wells in the case of 
clusters, instead of projected systems, these difficul- 
ties can be significantly alleviated. X-ray selection 
is less affected by projection effects, is done entirely 
by objective criteria and is subject to biases which are 
more easily quantifiable than those inherent in optical 
selection, particularly when an appropriate detection 
technique is used to identify extended X-ray sources 
and estimate their structural parameters. 

In this spirit, we have started a serendipitous 
search for galaxy clusters in deep (T > 20 ksec) 
ROSAT-PSPC pointed observations at high galac- 
tic latitude. The survey is being carried out using 
a new source detection technique for X-ray images 
based on the Wavelet Transform (Rosati, 1995). This 
approach has the unique capability to detect faint X- 
ray sources with greater sensitivity over a wide range 
in sizes and surface brightnesses when compared with 
standard detection methods; thus probing a larger 
range of cluster parameters. 

2. Observations and Data Analysis 

The data are 15 ROSAT-PSPC pointed observa- 
tions drawn largely from the public archive. Stellar 
targets or survey fields are selected with T > 20 Ksec 
and l&l > 20°. For AGN targets only cluster candi- 
dates more than 5' from the field center are retained. 

We have developed a fully automatic procedure for 
a homogeneous data reduction of the PSPC observa- 
tions. An image is extracted in the hard band [0.5-2.0 
Kev] with a small binning factor (8" pixel size) to ad- 



equately oversample the PSF. An energy-dependent 
exposure and vignetting correction map is applied. A 
Wavelet Detection Algorithm (WDA) is then used to 
localize and characterize sources. 

The Wavelet Transform (WT) has been exten- 
sively discussed in the literature (Grossman et al. 
1988, Bijaoui et al. 1991) and has been applied to 
a variety of astronomical problems (e.g, Slezak et 
al. 1993, Coupinot et al. 1992). The method de- 
veloped for the analysis of the X-ray images is de- 
scribed elsewhere (Rosati et al. 1993, Rosati 1995). 
Briefiy, it consists of an orthonormal decomposition 
of a signal into both space and scale through a con- 
volution with a family of functions called wavelets. 
The continuous WT of an image I(x,y) with respect 
to the analyzing wavelet tpaix, y) is the 3-dimensional 
set: W{x, y, a) = -^Hx, y) (g) t/'*(f ) |)) where a is the 
wavelet scale. The wavelet coefficient W(x, y, a) gives 
information about the signal at the location {x, y) for 
the scale a. In the Fourier space, the WT is inter- 
preted as a multi-scale filtering process. The signal is 
examined both in ordinary space, pixel by pixel, and 
in the frequency space, spatial scale by spatial scale. 
An inversion formula exists to recover the original im- 
age from its wavelet coefficients. 

A uniform sampling of the scale a in logarith- 
mic scale, flj- = (ao)', maximizes the orthogonality 
of the decomposition process. The wavelet used in 
this analysis is the difference of gaussians: ipair) = 

■^[e~~ — ^e~~]. This multi-scale analysis leads 
to an estimate of the morphological parameters (ex- 
tent and net counts) of the detected sources by fitting 
the wavelet coefficients with the WT of the model as- 
sumed for the sources. This fitting procedure is re- 
peated around each local maximum and the sources 
lying above the 99.95% confidence level (roughly 4(t) 
are retained. The final parameter estimate is obtained 
by weighting the structural parameters deduced from 
all the analyzed scales Oj-. The assumed source model 
is a gaussian isotropic profile which provides a good 
approximation of the PSF of the PSPC within an off- 
axis angle (6) of 15', hence only sources with 6 < 15' 
are retained. An anisotropic wavelet analysis (giving 
ellipticity and orientation parameters) is also possible 
and is used to provide a better evaluation of fiuxes and 
angular sizes of the cluster candidates. 

There are many advantages of the Wavelet Trans- 
form technique which make it particularly suitable for 
this project: i) the method is equally rigourous for 
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sources of different size; ii) no artificial parameters 
(such as detection/background box size) are needed; 
iii) morphological parameters can be estimated with- 
out knowledge of the background; iv) background gra- 
dients do not substantially affect detection efficiency 
and parameter estimation; v) sources are resolved 
down to the resolution limit and no blending is intro- 
duced; and vi) information on significant substructure 
within complex objects is preserved. 

In order to check the validity of our data analy- 
sis procedure, we show in fig.l the logA^-logS' of all 
the 601 point-like sources in the sample and note the 
excellent agreement with the results from a previous 
ROSAT deep survey (Hasinger et al. 1993, (H93)). 
It should be stressed that our survey reaches a fiux 
limit of 3 -10"^® erg cm~^ s~^ which is comparable to 
that of the II93 survey, while using shallower fields. 
In addition, no corrections for biases of the source 
detection procedure were required down to this fiux 
limit (we estimated these to be ^ 10% effect using 
simulations). This demonstrates the greater ability 
of the WDA over conventional detection techniques 
in characterizing sources in the low counts and high 
source density regime. 

For the conversion between count rate and fiux, a 
power law spectrum with energy index 1 for point-like 
sources and a thermal spectrum with T = 6 keV for 
extended sources have been assumed along with the 
galactic HI absorption appropriate for each field. 

The sky coverage has been computed as follows. 
For each field, the minimum detectable count rate in 
each radial bin of 1' width is calculated by integrating 
the average background over the average 95% power 
radius of the fitted profile and using the adopted de- 
tection threshold. For extended sources, the sensitiv- 
ity depends on the intrinsic extent as well. The lim- 
iting sensitivities and respective solid angles of these 
annuli for each of the 15 fields are then summed to 
obtained the cumulative survey area as a function of 
limiting fiux and source extent (see fig. 2). Above a 
fiux of ~ 2-10"" erg cm'^ s'^, the entire FOV (2.95 
deg^) of the 15 pointings contributes to the survey 
area for point-like sources. Complete sky coverage 
is reached at higher fiuxes for sources of increasing 
extent. 

3. Cluster sample selection 

Discrimination of extended emission, i.e. sources 
with characteristic size larger than the PSF, is more 
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Fig. 1. — Cumulative number counts for all the X- 
ray sources of the sample (721 sources with 6 < 18'). 
Error bars (Icr) are from this survey, solid line is the 
1(7 band from II93 data (adapted from Comastri et 
al. 1995). 




Fig. 2. — Sky coverage as function of fiux limit and 
intrinsic source extent computed considering the cen- 
tral 15' of each field. Measured fiuxes and intrinsic 
FWHMs of the 10 clusters in the sample are plotted 
on the top. 
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efficient in the centrai area of the detector where the 
PSF is sharper. In order to minimize the effects 
of confusion, in addition to retaining only extended 
sources with 6 < 15', we only used pointings shorter 
than 80,000 seconds. 

The FWHMs of the detected sources are compared 
with the best fit PSF as a function of the off-axis 
angle, as deduced from the whole set of evaluated 
sources. In order to identify possible systematic ef- 
fects in estimated fiuxes and sizes and to assess a sig- 
nificance for the sources to be extended, each field 
is simulated 10 times. This is done by replacing all 
the detected sources in the real fields with simulated 
sources drawn from a gaussian spatial distribution de- 
rived from the detected source parameters. The WDA 
is then re-run on the simulated images and the dis- 
tribution of the FWHM residuals of the PSF best fit 
model is compared with the one observed in the real 
fields (see fig. 3). The latter is a convolution of sta- 
tistical uncertainties in the parameter estimate with 
the FWHM intrinsic spread due to different source 
spectra and Njy absorption. From the distribution 
of FWHM residuals, one can derive a 99.9% c.l. for 
a source to be extended. Given the total number of 
detected sources, one can thus expect at most one 
source to be erroneously classified as extended as a 
result of random fiuctuations. The final selection cri- 
teria for the cluster candidate sample are as follows: 
1) Source FWHM > PSF FWHM at 99.9% c.l.; 2) 
fx >l - 10"" erg cm-2 s'^ = 

X-ray emission on scales larger than ~ 3' was not 
considered and hence corrections for background fiuc- 
tuations (e.g. Snowden et al. 1994) were not required. 
However, there is very little sensitivity for detection 
of these very extended sources at these fiux levels as 
can be seen from fig. 2. 

The fiux cut fmin allows us to maximize the com- 
pletness of the fiux limited sample as well as min- 
imize the contamination from confusion effect. Us- 
ing simulations with the observed log #-log S down 
to fiuxes of 10~^^ erg cm~^ s""'^, we estimate that at 
most 0.5 sources in our sample could be due to un- 
resolved emission of two blended sources below fmin ■ 
This low level of contamination is a result of the WDA 
performing adaptive aperture photometry, hence the 
confusion effect only becomes significant at fiuxes be- 
low ~ 5 • 10~^^ erg cm~^ s""'^. The same feature al- 
lows the fiux to be integrated over the largest radius 
consistent with the background (which is very low for 
the case of the ROSAT-PSPC). For the clusters with 



measured redshift we have verified this to be at least 
1 Mpc or approximately 4-5 core radii; thus limiting 
the lost fiux to 10-20%, well within the fiux errors. 
For the remaining clusters the R magnitudes of the 
brightest galaxies indicate a redshift interval between 
0.2 and 1, thus providing again a physical radius of 
~1 Mpc. 

The selection criteria yield a sample of 10 clus- 
ter candidates in 15 processed fields. The bright- 
est candidate has a fiux of 1.8 • 10~^^ erg cm~^ s~^ 
(~400 counts). Inspection of POSS-E plates show 
that the optical counterparts of these cluster can- 
didates are too faint to be readily identified. Pre- 
liminary CCD imaging of 7 candidates has revealed 
that at least 6 show a visible enhancement in the pro- 
jected galaxy density at the peak of the X-ray emis- 
sion. In addition, two further candidates turned out 
to be groups detected previously in other surveys (one 
at z = 0.34 in the Lockman field (M.Schmidt, pri- 
vate comunication) and the other at z = 0.13 in the 
Pavo field (Griffiths et al. 1983) ) and a third is a 
serendipitous rediscovery of a cluster in the Couch 
et al. (1991) distant cluster sample (J1888.16CL 
at z = 0.563). There are three more cluster can- 
didates visible in fig. 3 which have measured fiuxes 
(0.7 < /j; < 1) • lO""'^'* erg cm~^ s"-'^ of which two 
show a significant overdensity of galaxies around the 
X-ray peak. These preliminary optical follow-up ob- 
servations indicate that our technique has a high suc- 
cess rate. We show in fig. 5 [plate 1] two of the optical 
images and overlaid X-ray contours. 

4. Results: Cluster Number Counts 

With this sample, we have measured the surface 
density of clusters at fiuxes about an order of mag- 
nitude fainter than the limits reached by the Ein- 
stein Medium Sensistivity Sensitivity Survey (EMSS) 
(Gioia et al. 1990, Henry et al. 1992 (HY92)). We 
have derived a surface density of 6.7^2 q clusters per 
square degree at a fiux limit of 1 • 10"^'* erg cm~^ 
sec""'^. The errors correspond to the 68% confidence 
level and are deduced from the convolution of the 
dominating poisson statistics of the counts, N(>S), 
with the error distribution in N(>S) produced by the 
uncertainties in the estimated fiuxes and sizes as de- 
termined via simulations. In fig. 4, we have plotted a 
compilation of observed number counts over 4 decades 
in fiux. The number counts from the EMSS sam- 
ple have been obtained by integrating the luminos- 
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Fig. 3.— Estimated FWHM and 90% errors for all 
601 sources vs off-axis angle. The stars indicate clus- 
ter candidates; the solid line is a polynomial fit to all 
data points clipped at 2a. The inset shows the distri- 
bution of FWHM residuals. The 3 sigma lines were 
determined from Monte-Carlo simulations. 
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Fig. 4. — Observed and theoretical cluster cumulative 
number counts compiled from various sources. 



ity functions (XLF), n(L) = K{z)L~"^^\ reported in 
HY92 over the luminosity range (10'*^^ — lO'*®) erg s~^ 
in the three published redshift bins for z > 0.14 and 
using the Piccinotti et al. (1982) XLF for z < 0.14, 
i.e. a = 2.15, K = 3.8 • W'^ Mpc'^ L^^ (con- 
verted to the Einstein band and in agreement with the 
EMSS XLF within the systematic errors as discussed 
in HY92) (Hq = 50 km s"! Mpc'^, qo = 0.5 have 
been adopted). The integration is performed for two 
possible core radii (0.2 and 0.3 Mpc) which defines the 
closed area in fig. 4 once the fiuxes are converted to the 
ROSAT hard band (assuming a T = 6 keV thermal 
spectrum). Point 2 is from an X-ray selected sample 
of clusters from the ROSAT All Sky Survey (RASS) 
in the northern sky whose optical identification is 95% 
complete (Burg et al. 1994a). For reference, the de- 
rived number count from the Piccinotti sample is also 
plotted (point 3). 

We should note that the surface density of clus- 
ters at low fiux levels in fig. 4 should be considered a 
lower limit since the cut imposed on the extent ex- 
cludes clusters from the sample which are unresolved 
by the PSPC. Assuming a standard King profile for 
the surface brightness of clusters with /3 < 0.8 and 
core radius = 250 kpc, one can show that the 50% 
power radius exceeds 30" up to z = 1 in a critical uni- 
verse with Hq > 50 Km s~^ Mpc~^. Therefore, high 
redshift groups could possibly be missed in our sam- 
ple. The presence of cooling fiows with a S excess of 
10-30% in the center does not decrease the ability of 
the WDA to detect and classify sources as extended. 
However, extremely bright AGN in the center could 
result in a misclassification. The completeness func- 
tion could also be affected if cluster parameters such 
as Tc significantly change over the mass and redshift 
ranges probed. 

Results from ROSAT medium sensitivity surveys, 
whose optical identifications are now ~ 90% com- 
plete, can be used to place limits on the possible 
incompleteness of our sample due to misclassifica- 
tion of unresolved sources. These surveys, which uti- 
lize detection algorithms more geared to the detec- 
tion of point-like sources, are finding cluster number 
counts generally lower than our survey by a factor 5 
at /j; = 2 • lO""'^'* erg cm~^ s""'^ (Georgantopoulos et 
al. 1995). Hence, we expect the incompletness of our 
sample due to unresolved clusters to be less than 20%. 
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5. Discussion 

In fig. 4, we also overplot theoretical predictions for 
cluster number counts according to two models by 
Evrard & Henry, 1991 (EH91). The two curves have 
been shifted in the [0.5-2.0]keV band and rescaled to 
allow for the different normalization used for the XLF 
in the lowest redshift shell by a factor Keh9i/ Kz=o = 
7.2/3.8 = 1.9. These curves represent a CDM model 
(spectral index of the primordial fluctuations, n = 
— 1) with an _L cx M^(l + z)^/^ luminosity-mass re- 
lation (CDMP3) and a model with more power on 
large scales (n = -2) and L cx M"/S(l + z)"/-* 
(XCDM). Both models incorporate the negative evo- 
lution of the XLF as observed in the EMSS sample 
(see also Cavaliere et al. 1993 for a physical interpre- 
tation of alternative models). Also shown in fig. 4 is 
a curve representing a no evolution model which has 
been calculated by integrating the local (z < 0.14) 
XLF in the form of a power law down to luminosi- 
ties of 2.5 • 10'*^ erg""'^ s"-'^ and up to z = 1.5 using 
a = 2.15 and K = 3.8 • W'^ Mpc'^L^^ (see above). 
If a flattening of the XLF of X-ray selected samples is 
present, as has been observed in the XLF of optically 
selected samples (e.g. Burg et al. 1994b), then the 
same no evolution model would have a factor of two 
fewer counts at the survey flux limit and thus closer 
to the measured counts. 

At present, small number statistics and the un- 
certainty in evaluating the completness of our sam- 
ple prevent us from ruling out any of these models. 
More rigorous constraints will be possible once the 
redshift distribution and the volume density of our 
entire sample are known and compared with the re- 
sults from complementary studies, such us the RASS, 
which probe the XLF on a wider range of luminosities 
and at lower redshifts. Our ultimate aim is to use 
our deep survey to study the faint end of the XLF 
at moderate-high redshifts thus addressing a crucial 
missing link in our current understanding of cluster 
evolution. 
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RIXOS project. RB, RG and PR acknowledge par- 
tial support from NASA grants NAG5-1538, NAGW- 
2508, NAG8-794 and RDC acknowledges support 
from NAGW3288. 
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